CAS-KITPC/ITP-005 



An Ideal Channel of Long Distance Entanglement in Spin Systems 

Xiao-Qiang Xi 1 ' 2 , T. Zhang 3 , R. H. Yue 3 - 4 , X. C. Xie 2 , W. M. Liu 2 
1 Department of Applied Mathematics and Physics, 
Xi'an Institute of Posts and Telecommunications, Xi'an 710061, China 
2 Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100080, China 
3 Institute of Modern Physics, Northwest University, Xi'an 710069, China and 
4 Department of Physics, Ningbo University, Ningbo 315211, China 



o 
o 

(N 

so: 
< 

(N 

S : 
^— > ■ 

c ■ 
a ■ 

=5 

cr 



> 
m 
oo 

cn 

od 
O 

o 



X 



We propose a scheme for using spin chain to realize an ideal channel of long distance entanglement. 
The results show that there has different entanglement in different Hilbert subspace, the anisotropic 
parameter A will frustrate the entanglement and the magnetic field affect the entanglement through 
changing the ground state, the boundary entanglement C\n has the simplest expression in the 
simplest subspace and it only depend on the first item of the ground state, that item can be 
increased when a local magnetic field is introduced. Our propose can be handled easily because it 
only needs a uniform XX open chain initialized in the simplest Hilbert subspace and a bulk magnetic 
field that absent for the boundary qubits. 



PACS numbers: 03.65.Ud, 05.50+q, 75.10.Jm 

Introduction. —Quantum entanglement has played an 
important role in modern physics, for example, it is used 
to test some fundamental questions of the quantum me- 
chanics 1] and to act a central role in quantum informa- 
tion processing, like teleportation |2j, super-dense cord- 
ing @ j quantum computational speed-ups 0, § , quan- 
tum cryptographic protocols @, 0], one-way quantum 
computation [8[ and so on. It is also widely used in 
sensitive interferometric measurements [t| [lCj, 11, Hp 
and studying strongly correlated quantum systems [13j |. 
Especially, the ground-state entanglement can be re- 
lated to quantum phase transition [14] . Mott insulator- 
superfluid transition and quantum magnet-paramagnet 
transition. All of the applications of entanglement are 
closely dependent on how to produce it. Now, there are 
ma ny p hysical systems suggested to realize entanglement 
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tween the next-nearest-neighbor qubits in the transverse 
Ising model, but the maximal value of entanglement is 
about 4.3 x 10 -3 14[. A scheme is proposed for using 



Spin chain is a nature candidate for producing pair- 
wise entanglement. It has been used to quantum in- 

i. For the 



formation processing [21], 122, |23|, [2J, |2 



a five-qubit open chain with magnetic and system im- 
purities to realize a boundary entanglement with max- 
imal value of \ [1H. Another scheme is proposed for 
realizing an ideal boundary entanglement in four-qubit 
open chain with symmetry interaction [32| . Venuti et al. 
[331 ] propose a scheme with a "strong-weak-strong-weak- 
. . . " nearest interaction (dimerization) and a uniform 
next nearest interaction to realize a long-distance entan- 
glement in spin system and further propose it to qubit 
teleportation and transfer. The starting point of those 
works is analyzing the affect of the interaction to the en- 
tanglement, their different lies in their analytical method. 

The aim of this paper is to study how the interac- 
tion affect the non- nearest pairwise entanglement, espe- 
cially the boundary entanglement in the open chain. We 
will analyze the corresponding results then use them as 
guidance to construct a practicable ideal channel of long- 
distance entanglement. The Hamiltonian of Heisenberg 
XXZ open chain with impurity is 



pairwise entanglement, most of works focused on the en- 
tanglement between the nearest pair, which is basic and 
will help us to deeply understand the character of en- 
tanglement, and looked the ways to control and maxi- 
mize the entanglement, the effectively controlled factors 
include temperature, interchange coupling 23, |27|, mag- 
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netic field and system impurity [28|, [2 
entanglement can be realized for the nearest pair, while 
the nearest pairwise entanglement is not enough for the 
practical applications because its short distance. Gener- 
ally speaking, in solid system the entanglement between 
a pair decreases quickly as the increase of distance, so 
the works about non-nearest pairwise entanglement is 
scarce comparing with that of nearest, every non-nearest 
ideal entanglement will take great contribution to appli- 
cation of entanglement. There exist entanglement be- 



where Ji is the interaction between the i-th and (z+l)-th 
qubits, A is the anisotropic parameter, Bi is the magnetic 
3fJ. The ideal field, a ± = \(o x ±io y ), cr x ,a y , a z are the Pauli matrices. 



The anal ytic al tool in this paper is the concurrence 
theory 3J, [35|, Cij denotes the pairwise entanglement 
between the z-th and j-th qubits, which ranges from to 
1 is monotonous to the entanglement. 

Three-qubit Heisenberg XXZ open Chain in the uni- 
form Magnetic Field. —When the interaction Ji and 
the magnetic field Bi are uniform, the eigenvalues of 
the system are E , 7 = A =p 3B, Eia = ±B, E 2 ^ = 



-^±B,E 3 , 6 = ~^-±B, where A± = A±V8J 2 + A 2 
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and the corresponding eigenvectors are 



are uniform the eigenvectors can be constructed as 



IV-o > 
IV-i > 

\i>m > 
|V>4 > 

IV>„ > 

where C21 

C22 = c 52 = 

C32 = C 6 2 = 



|000 >, IV7 >= |m >, 
(-|001 > +|100 >)/v / 2, 

Cmi(|001 > +|100 >) + c m2 |010 > (m = 2,3), 
(-|011 > +|110 >)/v / 2, 

c„i(|011 > +|110 >) + c n2 |101 > (n = 5,6), (2) 



= C51 = 



^/2A 2 _ + 16J 2 



, C 2 2 = C 52 



4J 

A /2A T +T6J2 
4J 



C31 = C 6 i = 



4J 

^/2A 2 +16.7 2 ' 

and 



^/2A 2 + 16.7 2 



^A^ + ie./ 2 ' 

Using the concurrence theory, the entanglement be- 
tween the boundary qubits C13 can be obtained. The 
expression of Ciz(B, T, J, A) is tedious, the numerical re- 
sults show that the anisotropic parameter will frustrate 
the entanglement, and there exist a critical magnetic 
held for C13, this phenomena comes from the change 
of ground state as the magnetic field variety. When 
3A+V8J 2 +A 2 < B, the ground state is |V >> C i3 = °5 
when < B < 3A+V8 / 2+A2 , the ground state is |V> 2 >, 

(A-V8J 2 +A 2 ) 2 , , 

= — i — ' voy more clearly we 
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8J 2 + (A-V8.7 2 +A 2 ) 2 ' 

plotted the phase diagram of the ground state and the 
concurrence C\z{B) for some certain A: 
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Figure 1. (a) The phase diagram of the ground state, the shading 
is \ip 2 > and the blank is \ip >; (b) the concurrence Ci3(B) for 
different A. 

When the ground state is \ij) 2 >, C13 decreases as the 
increase of A, when A = 0, C13 = |, when A = 1, C13 = 
|, when A ^> J, C13 = 0. i.e. the magnetic field can be 
used to switch "on" and "off" the entanglement. From 
(b) of Figure 1. we see that the interval of B (in which 
the entanglement exist) increases as the increase of A, 
that is to say the large interval of entanglement existing 
is based on the decrease of entanglement. When B = 
0, the ground states are \tjj 2 > an d 1^5 > (duplicate 
degeneracy), C13 = max{2c2 1 — ^2,0} = 0. 

Four-qubit Heisenberg XXZ Open Chain in the uniform 
Magnetic Field. —As three qubit case, when Ji and Bi 



\if> > = |0000 >, IV15 >= |llll >, 
|V m > = c ml (|0001 > +e ia '» 1 |1000 >) 

+c m2 (|0010 > +e ia ™ 2 |0100 >) (m = 1, 2, 3, 4), 
|Vn > = c n i(|1110 > +e lQ "i|0111 >) 

+c„ 2 (|1011 > +e ia " 2 1 1000 >) (n = 5, 6, 7, 8), 
> = c fc i(|0011 > +e lQ "|1100 >) 
+c fe2 (|0101 > +e ta ^ 2 |1010 >) 
+c fe3 (|1001 > +e ia ^|0110 >) 

(k = 9, 10, 11, 12, 13, 14), (3) 

where the parameters c^, (0,7r) are determined by 
H\ip >= E\tp > and the normalization condition. 

The results of four-qubit case is similar as that of 
three-qubit, the magnetic field change the entanglement 
through changing the ground state. The entanglement in 
the ground state is the maximal. We picked out the pos- 
sible ground states and their eigenvalues: \ip > (E = 
~ 4_B), I Vi > (the state in \tp m > with minimal eigen- 
value Ei = -\{AB + J + V5J 2 + 2JA + A 2 )), \ip 9 > (the 
state in \ip k > with minimal eigenvalue Eg, which can be 
obtained for concrete A and J) and plotted the phase 
diagram of the ground state and the concurrence Cu{B) 
for some certain A in Figure 2. 
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Figure 2. (a) The phase diagram of the ground state; (b) the 
concurrence C\a{B) for different A. 

The maximal value of C14 for different A and ground 
state are shown in Table 1. 

Table 1. The entanglement C14 in the ground state, "GS(E)" 
instead of the ground state (eigenvalue). 



A 


B 


GS(E) 


Cl4 max 





< B < Vr >- 1 






0.0472 




V5— 1 ^ p ^ Vo+l 
4 ~- a ~- 4 




IV>i > 


0.2764 




75+1 < g 




\iPo > 





0.5 


< B < 0.48 


IV>9 


> (-2.712) 







0.48 < B < 1.25 




IV>1 > 


0.2 




1.25 < B 




IV-o > 





1 


< B < 0.66 


IV>9 


> (-3.232) 







0.66 < B < 1.70 




l^i > 


0.1464 




1.70 < B 




IV-o > 





2 


< B < 1.04 


IV>9 


> (-4.372) 


0.0149 




1.04 < B < 2.65 




IV>i > 


0.084 




2.65 < B 




IV>o > 
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When the ground state is >, C14 = 2c\ x will decrease 
as the increase of A. When the ground state is \ip g >, C14 
is complex as the variety of A: C14 = 0.0472 if A = 0; 
C14 = 0if0<A<l; C14 increases as the increase of A 
if 1 < A < 7, C14 decreases as the increase of A if 7 < A, 
the maximal value is about 0.063. 

Heisenberg XX Open Chain with system impurity as 
an "ideal" entanglement channel. —From the conclusion 
of three- and four-qubit cases, we see that the boundary 
entanglement obtain its maximal value at A = then de- 
crease as the increase of A. The uniform magnetic field 
can change the ground state and the corresponding en- 
tanglement. So in this section we no longer consider the 
anisotropy parameter A and the magnetic field Bi , XXZ 
model degenerate into XX model and we also suppose 
that all the degeneracy are eliminated. 

For convenience, we use total spin to sign the ground 
state in different Hilbert subspace, if the ground state 
has k spin up and N — k spin down, then total spin 
S T = - I = N ~2 k ■ For example, the total spin of 

IV'O >; \i>m >. IV'fe >, \i>n > and 1^15 > in Ec f ( 3 ) is 

2, 1, 0, —1 and —2 respectively. Because the symmetry of 
spin chain, we only consider the ground state with total 



„_;„ c _ iV iV -i N 

spin br — — , — -1, - — 

-iJV 



iV-2[#l 



ing dimension are C^ 
S T - 2 

in the other subspace, St 

N is the simplest, while St 



, the correspond- 
f— 1 

C}^ respectively. 



^ is a trivial subspace with dimension C$ = 1, 

1 with the minimal 



N_ 

2 



dimension C N 
is the most complex. The entanglement of the system is 
bounded by the ground state entanglement of the sub- 
space. 

Without A and B, the only interaction left in Eq. (1) 
is the exchange hopping J i; i.e. H = X^i* Ji{ a t a i+i + 
a i+i a T)- The case 01 N=4 is discussed thoroughly in 
Rcf. [? ]. Here wc study it in the Hilbert subspace. Let 
J± = J3 = 1, J2 = J, for the ground state with St = 1 
(i.e. \ip m > in Eq. (3), one spin up), Ci 4 (l) = 2c\ x — 

2(J-VI+J 2 ") 2 _ ./ 2 -./v^+j7 +2 
2(.7-yi+7 2 ") 2 +2(2) 2 
J 



. Cl4(l), 



when 



.7 2 -./V4+J 2 +4' ^^y'max 2 

0; for the ground state with St = (i.e. \ip k > 



in Eq. (3), two spin up), Ci4(0) = 2(|2cgiCg2| 

1 if J is large enough. 
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For N = 5, let J x = J 4 = 1, J 2 = J 3 = J. Calcu- 
lations show that Ci5(|) = 2 +\j* wu ^ S e t its maximal 
i when J — > 0, the analytical result of Ci5(|) can not 
be obtained, the numerical results show that Ci5(|) — > 1 
when J — > oo. 

When N=6, let J a = J 5 = 1, J 2 = J 3 = J 4 = J- The 
analytical result of C\§ at the ground state with St = 
2,1,0 can not be figured out. Numerical results show 
that C w (2) max —> \ when J -> 0; Ci 6 (l) TOOX -» 0.055 
when J w 2; Ci 6 (0) mm — 0.8 when J — > oo. If one let 
Ji = J5 = 1, J2 = J4 = J arid J3 = J 2 , Ci6(0) maa; = 1 
when J is large. 

Using this method, one can obtain the boundary en- 



tanglement of N-qubit in the case of a matrix with order 

rjv^fL "Vim. can t> e manufactured. In fact only N < 15 
I "2" J • ' ~ I "2" \ ) ■ 

(when N = 15, the matrix order is 429, great less than 
2 15 = 32768) case can be calculated. Although the ideal 
entanglement spin channel with interaction Ji = Jn-% = 
J 1 ^ 1 can be realized, this propose has great disadvantage 
when N is large, because realizing a very strong interac- 
tion is very difficult. This channel is practicable only for 
a short distance entanglement. 

A really ideal entanglement channel —A really ideal 
entanglement channel must be enough long and easy to 
manipulate. For the system with Eq. (1), the affect of 
A, Ji and the uniform magnetic field to the entanglement 
has been discussed in the above sections. They have no 
direct contribution to realize ideal entanglement channel. 
In the Hilbert subspace of N-qubit Heisenberg XX open, 
the simplest subspace is St = ~ 1 except for the trivial 
cast St — ^r, such a subspace can be really manipulated 
when N is enough large, so we choose such a subspace 
as a candidate for studying the boundary entanglement 
(Cin — 2c 2 nl ), the most important thing is to distinguish 
the component c m i(|l > +e Mml |iV >) (m = 1, 2, . . . , N) 
with the others, a good idea is to introduce the nonuni- 
form magnetic field, in which the simplest case is to ab- 
sent the magnetic field for the boundary qubits. So the 
Hamiltonian of the system can be written as 



N-l 



JV-l 



h = jJ2 k + ^+i + ^vn +bJ2 ( 4 ) 



The initial state of the system is prepared as the 
first qubit to be spin up and the others spin down. 
Since the Hamiltonian commutes with the total spin 
component along the z direction, the relevant Hilbert 
subspace must be spanned by the states \j >= 

|Oi,0 2 ,...,Oj_i,l :7 -,0 :H -i,...,Oiv > with j = 1,...,N. 
So the eigenvectors of the system can be written as 



\1>m> = J2 C "v(\i >+etam3 \ N + 1 -j > ')' N = 2k ( 5 ) 
J = l 

k 

\^m> = ^2c mj (\j>+e ia ^\N + l-j>) 



+c m , k +i\k + l >,N = 2k + l 



(6) 



where c mj , a mj ( or ir), j = 1, 2, . . . , k, m = 1, 2, . . . , N, 
are the parameters determined by H\ip >= E\ip > and 
the normalization condition. Because N — 2k is simpler 
than TV = 2k + 1, so we only consider even N case in the 
following discussion. 

The maximal entanglement is determined by the 
ground state, suppose \ip 1 > is the ground state, then 
C\n — 2c\ 1 will approximate to 1 if J2i=2 c ii ^ c ii- 

For N — 2k case, one need to solve two k x k matrix, 



4 



they are 



( x x J 

J X J 

J x 



M+ = 












y o o o 



x J 
J x J 
J x±J J 



(7) 



kxk 



with xi = -(2k - 2)B, x = -(2k- A)B and find its 

eigenvalues and eigenvectors. Only k=2 case can be 

calculated exactlv Cm - (2B- J+y-4i? 2 -4B./+5J 2 ) 2 
calculated exactly G 14 - (2i3 _,7 + v4B 2 -4i3j+5J 2 ) 2 +(2j) 2 > 

which is great than 0.99 if B/J > 5, C u = 0.2764 if 
5 = 0. When k > 3, one can only obtain the numeri- 
cal result of Cin , fortunately, the numerical results show 
that the parameters in the ground state \ip^ > satisfy 



\cij\ _ 2B 
\ c i,i + i\ 



2j- = ft, C\m = 2c^ 



_ /3 2fc (/3 2 -l) 



which is 



/3 2 (/3 2fc -l) ' 

great than 0.99 if (3 > 10. Details can be seen from Fig- 
ure 3. 



are out of the magnetic field), the long distance ideal en- 
tanglement can be realized in the simplest Hilbert sub- 
space. Our scheme needs two conditions, they are a uni- 
form XX open with interaction J initialized in the sim- 
plest Hilbert subspace and a bulk magnetic field B absent 
for the boundary qubits. Under these conditions C\m will 
be great than 0.99 if B/J > 5 for any even N. A uniform 
interaction chain is easier to realize than a chain with a 
"strong- weak-strong- weak-- • ■— " nearest interaction and 
a uniform next nearest interaction, in this aspect our 



scheme is simpler than that of Venuti et aUs [33[. If our 
scheme can be realized in experiment, then teleportation 
in a solid system will become reality. 
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Figure 3. (a) The diagram of C14 with B/J; (b) The diagram of 
C 1N with 2B/J = (3. 

In theory, the length of ideal entanglement channel can 
be infinite when 2B/J is large enough. As the increase 
of N, the eigenvalue different between the ground state 
and excited states will be smaller and smaller. 

Conclusions —The results for three- and four-qubit 
cases can tell us the affect of the anisotropic parameter A 
and Bi to the entanglement: A will frustrate the bound- 
ary entanglement; the affect of Bi to entanglement is to 
eliminate the degeneracy and change the ground state. In 
fact we can use same method to calculate N(<15)-qubit 
cases, the present results are still valid. 

For the Heisenberg XX open chain with nonuniform 
symmetry interaction Ji = Jjy-j = J 1-1 and J is large 
enough, the ideal entanglement can be realized in the 
most complex subspace, this conclusion has more theo- 
retical meaning than its application's, while it is a good 
candidate if one needs a not too long distance entangle- 
ment. 

For the Heisenberg XX open chain with uniform inter- 
action J and a bulk magnetic field (the boundary qubits 



[1] J. S. Bell, Speakable and Unspeakable in quantum Me- 
chanices: Collected Papers on Quantum Philosophy 
(Cambridge 1987). 

[2] C. H. Bennett et ai, Phys. Rev. Lett. 70, 1895 (1993). 

[3] C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 
2881 (1992). 

[4] P. W. Shor, SIAM Journal on Computing 26, 1484 

(1997) . 

[5] L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 
[6] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 
[7] D. Deutsch et ai, Phys. Rev. Lett. 77, 2818 (1996). 
[8] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[9] A. N. Boto et at, Phys. Rev. Lett. 85, 2733 (2000). 
[10] J. P. Dowling, Phys. Rev. A 57, 4736 (1998). 
[11] R. Jozsa et al, Phys. Rev. Lett. 85, 2010 (2000). 
[12] J. J. Bollinger et al, Phys. Rev. A 54, R4649 (1996). 
[13] J. Preskill, J. Mod. Opt. 47, 127 (2000). 
[14] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 

032110 (2002). 
[15] D. Bouwmeester et ai, Nature 390, 575 (1997). 
[16] Q. A. Turchette et al, Phys. Rev. Lett. 81, 1525 (1998). 
[17] A. Rauschenbeutel et ai, Science. 288, 2024 (2000). 
[18] F. Yamaguchi and Y. Yamamoto, Appl. Phys. A 68, 1 
(1999). 

[19] Y. Makhlin, G. Schon and A. Shnirman, Nature 398, 305 
(1999). 

[20] A. S0rensen et al, Nature 409, 63 (2001). 

[21] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 

(1998) . 

[22] A. Imamoglu et al, Phys. Rev. Lett. 83, 4204 (1999). 

[23] X. G. Wang, Phys. Rev. A 64, 012313 (2001). 

[24] S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 

[25] V. Subrahmanyam, Phys. Rev. A 69, 034304 (2004). 

[26] T. J. Osborne and N. Linden, Phys. Rev. A 69, 052315 



5 



quant-ph/0009060 



(2004). 
[27] M. C. Arnesen et al, 
[28] X. G. Wang, H. C. Fu and A. I. Solomon, J. Phys. A 34, 

11307 (2001). 
[29] X. Q. Xi et al, Phys. Lett. A 297, 291 (2002). 
[30] A. Abliz et al, Phys. Rev. A 74, 052105 (2006). 
[31] X. Q. Xi et al, Acta. Phys. Sin. 55 2036 (2006); 



[32] X. Q. Xi and W. M. Liu, Chin. Phys. 16, 1858 (2007). 
[33] L. C. Venuti, C. D. E. Boschi and M. Roncaglia, Phys. 

Rev. Lett. 96, 247206 (2006); ibid 99, 060401 (2007). 
[34] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[35] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 

(1997). 



